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ABSTRACT 


A fundamental  problem  In  plane  elasticity  Is  studied: 
to  determine  for  a given  opening  the  amount  and  distribution 
of  reinforcing  material  required  to  satisfy  certain  specified 
stress  field  criteria.  The  relationship  between  the  local 
stress  field  around  the  opening  and  the  reinforcement  distri- 
bution is  presented  explicitly  such  that  either  one  of  them 
can  assume  the  role  of  an  Independent  parameter.  The  problem 
Is  set  up  as  Muskhellshvll 1 's  first  boundary  value  problem  In 
which  a function  Is  derived  to  represent  the  loads  due  to 
reinforcement.  A solution  is  worked  out  for  a general  mapping 
function  which  can  describe  almost  any  opening  of  practical 
interest.  Numerical  results  and  other  implications  for  a 
reinforced  circular  opening  are  discussed  In  greater  detail. 

INTRODUCTION 

The  use  of  reinforcements  at  the  boundaries  of  openings  of  various 
shapes  is  common  practice  in  ship  and  aircraft  construction.  Beginning 
with  the  Kirsch  solution  of  a circular  opening  in  an  infinite  plate  under 
uniform  tension,  there  are  in  the  literature  a great  many  papers  which 
deal  with  the  boundary  value  problem  of  the  reinforcement  of  openings  in 
plates.  With  a few  notable  exceptions  [1,2]  most  researchers  formulate 
these  investigations  as  a stress/strain  boundary  value  problem,  i.e.,  the 
stresses  and/or  strains  are  specified  (constrained)  at  infinity  as  well  as 
at  hole  boundaries  and  the  resulting  stress  field  is  then  computed.  The 
approach  is  therefore  one  of  trial -and-error  in  which  several  designs  of 
reinforced  openings  must  be  studied  before  acceptable  stress  levels  are 
found. 

Mansfield  [1]  attempted  the  inverse  approach  of  specifying  the 
complete  stress  field  in  the  plate  and  then  seeking  a compatible  opening 
with  variable  reinforcement.  His  stress  criteria  required  the  original 
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stress  field  to  be  unperturbed  by  the  presence  of  reinforced  holes. 

Mansfield,  therefore,  used  the  phrase  “neutral  holes'*.  The  results  he 

produced  were  rather  limited  in  application,  as  some  of  the  reinforcements 

required  to  render  the  holes  neutral  were  not  very  practicable.  Recently, 

Bjorkman  and  Richards  [2]  studied  the  inverse  plane  elasticity  problem  of 

an  infinite  plate  containing  a hole  wnich  leaves  the  first  stress 

invariant  (o  + o ) of  the  original  stress  field  unchanged.  They  called 
X y 

such  holes  "harmonic  holes";  the  phrase  is  justified  in  that  o *a  is  a 

X y 

harmonic  function  which  remains  unaffected  by  the  presence  of  such  a hole. 

They  have  used  the  usual  Muskhel ishvi 1 i formalism  and  have  further  shown  ' 

that  an  ellipse  is  a possible  shape  for  a harmonic  hole  where  the  ratio  of  , 

major  to  minor  axes  bears  a direct  relationship  to  the  ratio  of  the  normal 

stresses  at  infinity  in  two  directions.  While  the  treatment  of  tht  problem 

by  Bjorkman  and  Richards  [2]  is  quite  elegant,  the  actual  application  is  ! 

again  somewhat  limited  due  to  the  shape  (elliptical)  and  unreinforced 

nature  of  the  holes. 

The  approach  followed  in  the  present  paper  builds  partly  upon  an 
earlier  paper  by  the  author  [3]  and  partly  upon  the  "harmonic  hole" 
condition.  The  problem  can  thus  be  formulated:  determine  the  distribution 
of  reinforcement  around  a specified  opening  in  an  infinite  plate  under  bi- 
axial loading  such  that  the  sum  of  the  principal  stresses  (or  the  first 
invariant  of  the  stress  tensor)  remains  unchanged. 

I 

BOUNDARY  CONDITIONS 

Figure  1 shows  an  opening  of  general  shape  which  is  reinforced  with  a 
variable  amount  of  reinforcement.  The  forces  acting  on  a boundary  element, 
ds,  of  the  reinforcement  are  shown  in  Figures  2a  and  2b.  If  this  boundary 
is  to  be  in  elastic  equilibrium,  the  following  equations  must  hold  true 
(in  curvilinear  coordinates,  Figure  2a): 


Figure  1 - A Reinforced  Hole  in  a Plate 


(a)  Curvilitwar  Coordinates 


P tin  \1/ 


(b)  Cartesian  Coordinates 


Figure  2 - Forces  on  an  Element  of  Reinforcement 


and  (In  Cartesian  coordinates.  Figure  2b): 

d(P  sin  Hi)  » t(o  dx  - T dy) 

J ^ 

d(P  cos  <<)  « 

These  equations  assume  that  the  reinforcing  member  does  not  resist 
bending  and  Is  distributed  thinly  along  the  hole  boundary.  If  the  Airy 
stress  function,  U,  Is  Introduced  In  Equations  (IB),  total  Integration 
and  summation  of  the  two  equations  gives 


ay 


where  c Is  an  arbitrary  complex  constant.  Let 


(2) 


z = z(c) 


(3) 


be  a mapping  function  which  transforms  the  actual  opening  In  the  z-plane 
to  a unit  circular  opening  In  the  c-plane.  Equation  (3)  can  be  used  to 
obtain 


(4) 


where  o Is  the  boundary  value  of  c at  the  opening.  Equations  (4)  and  (2) 
can  be  combined  to  give 


iu  + i ly.  = 

3x  3y 


(5) 


where,  due  to  equivalence  of  tangential  strain  at  the  opening  boundary  and 
In  the  reinforcement, 

^=^(og-voJ  (5A) 

The  boundary  condition  at  the  hole  boundary  can  thus  be  written  as 


(j)(o)  + ■ (!>'  (o)  + ii((o) 

z'  (a) 


o 


P 

t 


[UMl 

V 


+ c 


(6) 


where  ((i(o)  and  41(0)  are  boundary  values  of  the  functions  41(c)  and  4^(c)  at 
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I 
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4 


the  opening.  These  two  complex  functions  define  the  entire  stress  field 
by 


and 


a 

a 


o ♦ 
a 


2 i 


a6 


p'rrzT 


[2(c)  ♦'(c) 


♦•(c)) 


(7) 


where 


In  order  that  Equations  (7)  also  predict  the  stresses  at  infinity, 
the  complex  functions  ♦(;)  and  ♦(c)  assume  the  following  forms 


♦(c)  » S z(c)  ♦ ♦q(c) 
♦(c)  = 0 z(c)  + vo(c) 


(8) 


where  S and  D are  constants  and  the  functions  ♦q(c)  and  ♦q(c)  approach 
zero  at  large  values  of  c-  Thus 


« a' 

'o'‘>  • 7 

n=0  c 


(9) 


The  constants  S and  0 can  easily  be  computed,  using  Equations  (7),  to  be 

(10) 


^ 4 


D = . 

where  p and  q are  the  stresses  at  infinity  and  0 is  the  angle  between  the 
x-axis  and  p.  The  form  of  ♦(c)  and  ♦(c)  can,  thus,  be  effectively  used  to 
automatically  satisfy  the  boundary  condition  at  infinity  and  Equation  (6) 
would  then  be  the  only  remaining  condition  to  be  satisfied  for  obtaining  a 
solution  to  the  problem.  It  should  be  noted  that  Equation  (6)  is 
considerably  simpler  than  a comparable  equation  developed  earlier  by  the 
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author  [3]  for  solving  the  direct  problems  of  elasticity  involving 
reinforced  hole  boundaries. 


METHOD  OF  SOLUTION 

Once  the  boundary  condition  as  given  by  Equation  (6)  is  developed, 
its  solution  is  merely  routine  and  may  or  may  not  be  complicated  depending 
on  the  complexity  of  the  mapping  function.  In  this  case,  however,  this 
equation  is  further  simplified  by  specifying  a field  condition  which 
requires 


p + q 


(11) 


Equation  (11)  reduces  to  Oq(c)  = 0,*  as  shown  by  Bjorkman  and 
Richards  [2].  Thus  Equation  (6)  becomes 


2S  z(o)  * 0 z(o)  + i<)q(o)  ° q ^ ^ z'  ('P 

If  the  opening  is  assumed  to  be  a snwoth  closed  curve,  it  is  possible 
to  represent 

I (co"  + ^)  (13) 

Since  (o)  z'  (a)'  is  a real  quantity,  c^  and  d^  have  to  be  accordingly 
adjusted.  Equations  (12)  and  (13)  can  be  combined  to  obtain 


2S  z(o)  z‘  (o)  + D z(o)  z'(o)  + iJ»Q(a)  z'(o)  = 

E (c  + ^)  + c (14) 

n=0  0 

which  is  now  the  boundary  condition  at  the  reinforced  opening  of  any 
general  shape.  It  is  easy  to  show  that,  with  Muskhelishvili 's  technique. 


* Care  must  be  taken  since  Bjorkman  and  Richards  [2]  use  different  notation. 


L. 


6 


rt-|-0 


Equation  (14)  can  conveniently  be  solved  for  a general  mapping  function  of 
the  type 

2(c)  » R(c  + £ ^ ) (15) 

m=l  c 

which  can  represent  almost  any  opening  of  practical  interest.  However, 
to  focus  attention  on  a particular  opening  and  to  study  the  related 
implications,  the  problem  of  a circular  opening  will  now  be  discussed  in 
greater  detail . 


CIRCULAR  OPENING 


Let  a plane  plate  under  biaxial  uniform  tensions  p and  q in  the  x- 
and  y-directions,  respectively  (see  Equations  (10)),  contain  a circular 
opening  of  radius  R.  The  mapping  function 


z = ziz)  « Rt  ; c = re 
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(16) 


which  describes  a circular  geometry,  can  be  substituted  in  the  boundary 
condition  (Equation  (14))  to  obtain 

d 


2S  R2o  + — + R ! b o"  = Me  o"'*'’ 


n=0 


n=0 


n-1 


) +cR 


(17) 


Now  Equation  (17)  and  its  conjugate  can  be  multiplied  by 

ZlTl  o-c 

and  integrated  around  the  unit  circle  to  obtain 


Cn  = dn;  n M 
2Cq  = b^R  + 2S  R2 

"l  = ° 

C2  = 0 R2 

c^  = 0;  n > 3 


c = 0 


b2  = 0 

bj  = DR 

b^  = 0;  n > 4 


(18) 


which  define  the  solution  to  the  problem.  The  functions  (|>(c).  and 

can  now  be  given  by 
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1 


♦ (c)  = SRc 

♦ (c)  = DRc  + ^ ^ 

C 

^ = b,  + 2SR  + DR(o2  + 4-) 


Equations  (5A),  (7),  and  (19)  can  be  combined  to  obtain 


o = 2S  + -V  - (0  - % cos  2e 
r Rr2  r**' 


»e  ' - 572-  * 


(19) 


(20) 


re 


= (D  + ^)  sin  2e 


and,  since  at  the  hole  boundary  r Is  unity, 

^ ^ (p+q)  +2  - 2(p-q)cos  2e 

_ = 


(21) 


(1-v)(p+q)  - 2(l+v)  + 2(l+v)(p-q)cos  2e 

Equations  (20)  and  (21)  clearly  show  a coupling  through  the  term 

^1 

Actually,  the  presence  of  this  term  Implies  that  several  pairs  of 
reinforcements  and  stress  distributions  are  possible  for  a reinforced 
circular  opening  which  do  not  alter  the  sum  of  principal  stresses  In  the 
rest  of  the  Infinite  plate  under  uniform  biaxial  tension  at  Infinity.  The 

^1 

possible  values  of  p,  q,  and  , which  do  not  permit  the  actual  rein- 
forcement distribution  A/Rt  to  become  negative,  are  given  by  the  following 
Inequalities: 

3 + 2 

'’  ‘^qR  q 3 3qR 


3+v  p 1+  V ^ ^ ^ H3V  ^ 2 l'*'v  ^1 
l+3v  l+3v  qR  q 3+v  3+v  qR 


(22) 


These  equations.  If  plotted  on  graph  paper,  will  delimit  an  area  which 
will  define  several  permissible  pairs  of  values  for  p/q  and  bi/qR.  It  Is 
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obvious  that,  for  p=q.  Equation  (21)  reduces  to  the  Mansfield  solution  [1]. 
At  this  point  it  should  be  noted  that  right  from  the  beginning  the  problem 
did  not  call  for  a unique  solution.  The  constraint  of  unaltered  sum  of 
principal  stresses  does  not  characterize  the  stress  field  as  unique.  As  a 
result  a number  of  possible  solutions  exist  which  satisfy  the  prescribed 
conditions.  Any  additional  condition,  such  as  regarding  the  stresses  at 
any  particular  location,  will  immediately  make  the  solution  unique.  The 
general  procedure,  as  described,  allows  for  some  trade-off  studies 
between  the  reinforcement  distributions  and  the  stresses.  Of  course,  as 
pointed  out  by  Bjorkman  and  Richards  [2],  other  field  conditions  can  also 
be  prescribed  as  the  beginning  point. 

To  get  an  actual  feeling  for  the  solution,  let 

n b, 

^ = 1.5  and  4 = -0.5 
q qR 

Then 

^ _ 3-2  cos  2e 

Rt  (7-3v/  + 2(l+v)cos  26 

which  has  a maximum  value  of  at  e = it/2  and  a minimum  value  of 

l-v  9-v 

at  s=0.  This  distribution  obviously  requires  considerably  less 
reinforcing  material  than  the  Mansfield  case  of  all-around  tension  with 
p/q  = 1 and  a uniform  A/Rt  = 1/1-v. 

ADDITIONAL  EXAMPLES 


Circular  geometries  are  a special  case  in  many  physical  problems  and 
their  solution  is  relatively  simple.  Therefore,  it  appears  worthwhile  to 
investigate  some  noncircular  openings.  Let  an  infinite  plate  under  uni- 
form biaxial  tensions  p and  q along  x-  and  y-axis,  respectively,  contain 
a reinforced  opening  given  by 


(23) 


It  is  well -documented  that  a function  such  as  Equation  (23)  can  map 
a circular,  elliptical,  square,  or  rectangular  opening  in  the  z-plane  onto 
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a unit  circle  in  the  c-plane  depending  on  the  values  of  constants  a^ , a^, 
and  a^.  Furthermore,  if  the  orientation  of  the  opening  is  symmetric  with 
respect  to  the  axes,  then  the  constants  a^ , a^,  and  a^  are  real.  Now  with 
Equation  (23)  as  the  mapping  function,  the  boundary  condition  (14)  after 
usual  processing  yields 

c_  = d ; n ^ 0 
n n 

2Cq  = Rb,  + 2SR2(l-a2-3a2-5a2) 

Co  = 2SR2(a,-a,ao-3aoac)  + DR2 
2 11335  (24) 

c^  = 2SR2(a2-a^a5) 

Cj  = 2SR2  aj 

c^  = 0;  n > 7 

The  function  lii(c)  can  be  determined  by  integrating  the  conjugate  of 
Equation  (14).  In  fact,  the  coefficients  are  found  to  be: 

b2^  = 0;  n > 0 

Cp 

bj  = a^b^  X ^ 2SR(a^+3a^a2+5a2a5)  + DR(a|+2a2) 

c. 

bg  = a^b2+3a2b^+ -^+2SR(3a2+5a^ag)+4DR(a^a2+ag) 

Cf-  , 

by  = a^bg+3a2b2+5agb^+ -^+10SRag+3DR(2a^aj+a2) 

(25) 

bg  = a-jby  + 3a2bg  + Ba^bj  + SDRa^a^ 

bll  " a^bg + 3a2by + Sa^bg + 5DRag 

^13  ®1^11  ^^3^9  ^^5^7 

^2n+l  " ^1*^20-!  ■^^®3^2n-3^^^5^2n-5’ 


Now  Equations  (lA)  and  (5A)  and  the  field  condition  can  be  combined 
to  obtain  a useful  relation  that  will  bypass  the  use  of  i/)(c);  i.e.,  for 
computing  A/t: 
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where 


A 

t 


P ds 

T de 

t - ('*'’>  F 


(26) 


(z'z‘)^^^ 


- z'z" 

I _ I i ^ * 


Z Z' 


+ — 7*7" 

2o  2 ^ ^ 


is  the  radius  of  curvature  of  the  opening  boundary.  The  problem  is  thus 
solved.  In  fact,  a direct  substitution  of  the  required  quantities  in 
Equation  (26)  leads  to 


where 


and 


A , P(o)(z'z')^^^ 

^ 4S(z'z')2-  (l+v)P(o)F(o) 


P(a) 


3 

I 

n=0 


^2n 


(26A) 


F(o)  = z'z'  + 


The  mapping  function  z is  given  by  Equation  (23)  and  the  constants  c^ 
by  Equation  (24).  The  situation  is  again  similar  to  that  of  the  circular 
case.  A/t  is  known  within  a constant  b^  which  couples  it  with  the  stress 
field.  Again  a number  of  stress  fields  can  be  paired  with  appropriate 
reinforcements  which  will  keep  the  sum  of  principal  stresses  unaltered 
over  the  entire  field. 

Using  Equation  (23)  for  a mapping  function  and  assuming  no  rein- 
forcement, the  problem  can  be  reduced  to  that  of  reference  [2],  which  is 
to  be  expected.  In  contrast  to  reference  [2],  which  determines  the 
harmonic  opening  shape  to  suit  an  existing  stress  field,  this  paper 
proposes  a method  for  reinforcing  an  opening  of  any  desired  shape  in  such 
a way  that  the  reinforced  opening  behaves  as  a "harmonic  hole." 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

(1)  DTNSRDC  REPORTS,  A FORMAL  SERIES  PUBLISHING  INFORMATION  OF 
PERMANENT  TECHNICAL  VALUE,  DESIGNATED  BY  A SERIAL  REPORT  NUMBER 

(2)  DEPARTMENTAL  REPORTS,  A SEMIFORMAL  SERIES,  RECORDING  INFORMA 
TION  OF  A PRELIMINARY  OR  TEMPORARY  NATURE,  OR  OF  LIMITED  INTEREST  OR 
SIGNIFICANCE,  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION 

(3)  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  USUALLY  INTERNAL 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS,  NUMBERED  AS  TM  SERIES 
REPORTS;  NOT  FOR  GENERAL  DISTRIBUTION. 


